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Abstract 



D , We consider constraints on CP- violating phases in the Constrained Minimal Supersymmetric 
. 

Standard Model. We find that by combining cosmological limits on gaugino masses with 
experimental bounds on the neutron and electron electric dipole moments, we can constrain 
^ I the phase of the Higgs mixing mass to be |6'^| < vr/lO, independent of choices of the other 
mass parameters in the model. The other CP-violating phase 6a is essentially unconstrained. 



The difficulties associated with trying to constrain the vast parameter space of the general 
low-energy Minimal Supersymmetric Standard model are well known, and some simplifying 
assumptions are usually be made in order to get an experimental handle on the set of new 
parameters. A common ansatz is that of the Constrained Minimal Supersymmetric Standard 
Model, based on supergravity and grand unification (described in more detail below). The 
CMSSM is interesting because it is simple, predictive and naturally provides a stable dark 
matter candidate (an LSP bino-type neutralino) over most of its parameter space. In [Q, 
we showed that there was a strong correlation between the CP-violating phases in the MSSM, 
the cosmological LSP relic density and the neutron and electron electric dipole moments, 
when all of the parameters are set at the weak scale and hence are independent of any RGE 
evolution. Here we wish to consider the effect of the CP-violating phases present (though 
normally ignored) in the CMSSM, specifically the inducing of electric dipole moments for 
the neutron and the electron. We will find that by combining cosmological constraints on 
the mass of the LSP with experimental bounds on the neutron and electron EDM's, we can 
bound one of the two new phases in the CMSSM to lie within \9^\ < vr/lO. 

The MSSM contains a plethora of new parameters, which make empirically testing the 
model difficult. In addition to the supersymmetric Higgs mixing mass parameter /i, there are 
supersymmetry-breaking gaugino masses M,, sfermion mass^ parameters m| , trilinear scalar 

Ji 

parameters Ai, the Higgs scalar mixing mass^ Bfi, and the ratio of the Higgs vevs, tan/5. 
This large parameter space can be simplified at the unification scale by taking a common 
gaugino mass M, sfermion mass^ parameter mg and trilinear scalar parameter A, and this 
model is referred to as the CMSSM. Since in the CMSSM, the two Higgs mass^ parameters 
are set equal to the sfermion mass^ parameters at Mx, the other Higgs sector masses are 
determined by the requirement of correct electroweak symmetry breaking, once tan (3 and 
rritop are fixed. In principle, M, fi, Bfi, and A may be complex; however, not all of these 
phases are physical 0]. It is possible to rotate away the phase of the gaugino masses. And, 
by making Bfi real, we ensure that the vacuum expectation values of the Higgs fields are 
real. The CMSSM, then, is specified by three masses (mo,M, and Aq), two phases (6'^(Mx) 
and 9a{Mx)), tan/3 and rritop- 

Once the masses and phases are given at Mx, they can be RGE evolved to the electroweak 
scale. In practice, we use the one-loop RGEs for the masses and two-loop RGEs for the gauge 
and Yukawa couplings Q]. The structure of the equations for the gauge couplings, gaugino 
masses and the diagonal elements of the sfermion masses are such that they are entirely real. 
The evolutions of the Ai, however, are more complicated, as the Ai pick up both real and 
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imaginary contributions. For example, the evolution of At is given by 

^ = ir^{-f9lMs- 39I M,-f gl Ml + hlA, + 6/^^*) (1) 

As one can see, At receives real contributions qM proportional to the gaugino mass (whose 
coefficients Cj are different for each sfermion in a generation) and (in principle complex) 
contributions dih'jAf from the heavy generation (whose coefficients di differ for the ffist two 
and the third generations); the phases (and magnitudes) of the Ai must therefore be run 
separately. At one loop, the evolution equation for fi is given by 



d/i /i 



{-^gl - gl + hl + 3hl + 3h^) (2) 



dt 167r2 

and the phase of /i does not run. 

The aim of this paper is to combine cosmological bounds on the relic density of neutralinos 
with experimental bounds on the neutron and electron electric dipole moments in order to 
place limits on the new CP- violating phases and 6a- The relic cosmological density of 
neutralinos and the electric dipole moments are both strongly dependent on the sfermion 
masses. We take the general form of the sfermion mass^ matrix to be [S] 



(3) 



fMl + m) + cos 2(3{Ty - Qf sin^ ^vf)M| ruf rUfe'^f \ 

V mffUfe-'^f Ml + m'j + cos2f3Qfsm^9wMl) 

where M^pt) are the soft supersymmetry breaking sfermion mass which we have assumed 
are generation independent and generation diagonal and hence real. Due to our choice of 
phases, there is a non-trivial phase associated with the off-diagonal entries, which we denote 
by m^(m^e*'^^), of the sfermion mass^ matrix, and 

nife'^f = Rffi + A} = Rf\i^\e'^^ + \Af\e''^^f , (4) 

where rrif is the mass of the fermion / and Rf = cot (3 (tan/5) for weak isospin +1/2 (- 
1/2) fermions. We also define the sfermion mixing angle ^/ by the unitary matrix U which 
diagonalizes the sfermion mass^ matrix, 

/ cos 9f sin 9f e^^f \ 
U=i . ' . (5) 

V — sin 6f e ^'^f cos Of ) 

Previously , it has been shown that the presence of new CP- violating phases may have 
a significant effect on the relic density of bino-type neutralinos. The dominant channel for 
bino annihilation is into fermion anti-fermion pairs. However, this process exhibits p-wave 
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suppression, so that the zero temperature piece of the thermally averaged annihilation cross- 
section (which is relevant for the annihilation of cold binos) is suppressed by a factor of 
the final state fermion mass^. This significantly reduces the annihilation rate and increases 
the neutralino relic density. Mixing between left and right sfermions lifts this suppression 
to some extent by allowing an s-wave contribution to the annihilation cross-section which 
is proportional to the bino mass^, but the presence of complex phases in the off-diagonal 
components of the sfermion mass matrices dramatically enhances this effect. Explicitly, we 
compute the relic density by using the method described in ref. and expand {aVrei) in a 
Taylor expansion in powers of T/m^ 



{avrei) = a + b (T/m^) + O (^(T/mg)') 



The coefficients a and b are given by 



/ 

-3 + 



(6) 



(7) 



3m} 



where df and 6/ are computed from the expansion of the matrix element squared in powers 
of p, the incoming bino momentum, and Vj = {1 — m'j/rn~y^'^ is a factor from the phase 
space integrals. For M| ^ and rrif <^ m~ g]. 



a/ 



327r 



5^ ^ Yl Yl sin2(2^/) sin^ 7. 



0(m^.m^) (9) 

We see that if sfermion mixing is significant {6f is large) and the phase of the mixing term 
is large, the p-wave suppression is removed. 

The effect of CP- violating phases in a simple model where all the sfermion scalar mass pa- 
rameters were (for convenience) taken equal at the electroweak scale, along with the trilinear 
A parameters, was studied in [|2| . It was found that the cosmological upper bound on the bino 
mass in this model was increased from 250 GeV [0, |^ to 650 GeV by the presence of the new 
phases. This is a significantly stronger effect than the enhancement due to mixing alone 0. 
The enhancement found in this case can be partially traced to the assumption of equal scalar 
masses at the weak scale. In particular, sfermion mixing is sensitive to this assumption as 
can be seen from the magnitude of sin^(26'/); for niffrif <^ M| — M|. + 2Q/ cos 2/5 sin^6'vK^, 



Z5 



sm 



2 2 



{Ml -Ml + 2Qf cos2/3sin2^vFM|)2 



(10) 
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The assumption of equal masses translates into taking Ml = Mr and hence the enhancement 
in s\.v?{26f). Constraints from bounds on the neutron electric dipole moment, however, 
restrict the range of the relevant phases to < vr/lO and 7down < 7i"/6, so that the mg(max) 
is reduced to about 350 GeV. It should be noted however, that the constraints from the 
neutron electric dipole moment were based on the naive quark model and may in fact be 
overly restrictive when the strangeness content of the nucleon is taken into account |1C]. 

The above mass ansatz of equal masses at the weak scale is particularly simple. However, 
this pattern of low-energy sfermion masses contains potentially dangerous charge and colour- 
breaking minima in the scalar potential l^- To look at the scalar potential above the 
electroweak scale, one must RGE evolve the mass^ parameters up to the scale of interest. 
In the MSSM, the sfermion mass^ parameters tend to run negative at large scales due to 
sfermion couplings to gauginos, and for some initial low-energy values of the scalar and 
gaugino masses, CCB minima appear at high scales. These CCB minima may be avoided by 
taking the low-energy sfermion masses to be sufficiently high. However, these large sfermion 
masses lead to a neutralino relic density which is much too large, unless both the common 
A parameter is tuned so that the lighter of the two stops has a mass close to the neutralino 
mass and the neutralino is heavier than the top. In contrast, the CMSSM is free of these 
minima by construction (although other CCB minima may be present [0]). In this case the 
sfermion mass^ parameters are taken non-negative and equal at the unification scale Mx and 
are driven more positive by the gaugino couplings as they are run down to the electroweak 
scale. In a fully consistent model, one of the Higgs mass^ parameters runs negative at low 
scales to provide SU{2) x U{1) breaking. 

As we have said above, in addition to its simplicity, the CMSSM also enjoys the feature 
that in most of the parameter space, the lightest neutralino x? (and in fact the LSP) is mostly 
binofll], and for large values of the unified gaugino mass M, x5 is almost pure bino. We can 
thus compute the neutralino annihilation rate including only fermion anti-fermion final states 
(as in (P)), and the error made in not considering higgsino mixing will be insignificant for 
the relevant gaugino masses. In Figure 1, we show contours of constant neutralino relic 
density, in the M-uiq plane. Throughout this paper we use tan/? = 2.1 and mtop = 170 GeV. 
Three contours are shown, for Vt-h"^ of 0.25, 0.5, and 1. The two vertical contours are 
bino purity contours of 0.95 and 0.99; we see that for bino masses near their upper bounds 
for ^^h'^ < 0.25, where mg ^ 0AM ^ 160 GeV, the higgsino admixture in the LSP is 
extremely small. The shaded regions in Figure 1 are ruled out because they lead either to a 
chargino with a mass less than ~ 65 GeV |T^, a sfermion with a mass less than 74 GeV[|T5|, 
or a stop, chargino, or (in the bottom right region) stau as the LSP. The contours of ^zh"^ 



4 



were computed for the particular choices of Aq = 300 GeV,^^ = 0,6*^ = 0.8 vr; however, the 
contours vary httle as these parameters are changed, and the variation is aheady completely 
negligible at the upper range of the allowed bino mass for ^^h? = 0.25. There is a similar 
picture for 6^ = ti {ji ^ —f^)- 

In contrast to the model studied in 0, the upper bound on the bino mass in the CMSSM 
is essentially independent of the CP- violating phases and 9 a. Recall from (^) that the 
phases are only important if there is a significant amount of L-R sfermion mixing. In the 
CMSSM, when M| and M|, are run down from Mx, they are split by roughly 0.4M^ at 
the electroweak scale. From an examination of Eq. (10), one can see that the off-diagonal 



components of the sfermion mass matrices which contain powers of the fermion mass are much 
less than M| — when M is a few hundred GeV, therefore sfermion mixing is negligible, 
and the phases cannot lift the p-wave suppression. Additionally, since at the upper range of 
bino masses we are deep in the pure bino region, the composition of the lightest neutralino 
is independent of /i, and 6^ in particular. One can therefore set an upper bound for M from 
cosmological considerations which is quite independent of parameter choices. (If one takes 
tan/3 very large so that stau mixing is significant, then the neutron and electron electric 
dipole moments given below become very large as well. Constraints on the CP-violating 
phases from bounds on EDM's then ensure that effect of the phases on the neutralino relic 
density are negligible.) One finds that for ^^h"^ < 0.25, M must be less than about 400 GeV, 
which corresponds to a bino mass of roughly 160 GeV. 

We turn now to the calculation of the electric dipole moments of the neutron and the 
electron. As our calculation of the EDM's parallels our previous calculation 0, we refer the 
reader there for details. As in @, we will only include the three contributions coming from 
neutralino, chargino, and gluino exchange to the quark electric dipole moment. The necessary 
CP violation in these contributions comes from either 7j in the sfermion mass matrices or 6*^ 
in the neutralino and chargino mass matrices. Full expressions for the chargino, neutralino 
and gluino exchange contributions are found in [0. 

The contributions to the quark electric dipole moments from the individual gaugino ex- 
change diagrams fall as M is increased, because the squark masses^ receive large contributions 
proportional to during their RGE evolution from Mx to Mz- Roughly, 

m|^mg + 6M2 + 0(M|). (11) 

Thus even for large values of the CP violating phases, one can always turn off the quark 
electric dipole moment contributions to the neutron EDM by making M sufficiently large [|16|; 
however one must still satisfy the cosmological bounds discussed above. Experimental bounds 
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are \dn\ < 1-1 x 10 ^^e cm [|17| for the neutron electric dipole moment and \de\ < 1.9 x 



lO^^^e cm [|T^ for the electron EDM. Note also that the squark masses m~ are only weakly 
dependent on mo in the cosmologically allowed regions of Figure 1, and so the quark EDM's 
will also be independent of uiq. 

We have computed the neutron EDM in the CMSSM as a function of 6'^, Oa-, and M for 
fixed Aq, mo, and tan/5, using the naiVe quark model. In practice we find that the dominant 
contribution to the quark electric dipole moments in the CMSSM come from the chargino 
exchange diagrams, unless 6^ is extremely small (and 6^ ^ 6^). We can then find, as a 
function of 6'^ and 6a, the minimum value of M required to bring the quark electric dipole 
moment contributions to the neutron EDM below the experimental limits. In Figure 2, 
we plot contours of Mmin for Aq = 300 GeV and ttiq = 100 GeV. The light central region 
corresponds to Mmin < 200 GeV, and successive contours represent steps of 100 GeV. The 
darkest regions on the left and right sides of the figure lead to a stau as the LSP and 
correspond to the lower right shaded region in Figure 1; their positions do depend somewhat 
on mo, as can be seen from Figure 1. The EDM's were computed on a grid with spacings 
of tt/IO in 6a and tt/SO in 6'^, so contour features smaller than these dimensions are not 
significant. Recalling from Figure 1 that ^^h"^ < 0.25 requires M < 400 GeV, we find the 
constraint that 16*^1 < 7r/15 (the region with 6^ < comes from values of 6a near 37r/2). 
There is a similar allowed region near ^^^j = tt, corresponding to a sign change in /i. The 
dependence on 6a in Figure 2 is weak because it affects only the gluino exchange contribution, 
which is sub dominant, and even then 6a is only partially responsible for the relevant phases 
7up and 7down (see Eq. (|)). The figure is shifted slightly to one side of = because in 
one direction there is a cancellation between the chargino and gluino exchange pieces, and 
in the other direction the two contributions come in with like signs. These bounds are not 
dependent on mo and are only very weakly dependent on Aq for Aq < 1 TeV. 

In ref it was found that parameter choices which produced a sufficiently small neutron 
EDM also produced an electron EDM well below experimental bounds, so that it was never 
necessary to separately consider bounds from the electron EDM. In the case of the CMSSM, 
this is not always the case. The reason is that in the CMSSM, the squarks receive a large 
contribution to their masses from the RGE running and are consequently much heavier than 
the sleptons, in contrast to the case studied in ref p| where all the sfermions shared a common 
scalar mass^ parameter at Mz- The analysis for the electron EDM is slightly complicated by 
the fact that the slepton masses, and consequently the electron EDM, do depend somewhat 
on mo in the cosmologically allowed region, even at the upper limit of the neutralino mass, 
where M takes its largest allowed value. However, for values of mo below 100 GeV, this 
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dependence is not strong, and the minimum values of M we find change by less than 10% for 
values of mo < 100 GeV. In Figure 3 we show contours of constant Mmin for Aq = 300 GeV 
and rriQ = 100 GeV, with the requirement that \de\ < 1.9 x lO^^^ecm |]T^. Again, the darkest 
regions on the left and right sides of the plot lead to a stau as the LSP and correspond to 
the lower right shaded region in Figure 1. This time the slight dependence on 6*^ comes from 
the dependence of the neutralino exchange contribution (which in this case is significant) on 
the phase 6 a- We observe that we get comparable bounds to those from the neutron EDM 
in Figure 2, \6^\ < vr/lO. It is important to note however, that the bounds coming from the 
electron electric dipole moment are clearly not sensitive to the spin structure of the nucleon 

In summary, we have combined cosmological bounds on gaugino masses with experi- 
mental bounds on the neutron and electron electric dipole moments to constrain the new 
C7P-violating phases in the Constrained Minimal Supersymmetric Standard Model. We find 
that in contrast to models studied previously, the phases do not affect the cosmological limits 
on the mass of an LSP bino-type neutralino. While there is no bound on the phase 6^ of 
the unified trilinear scalar mass parameter A, the phase of the supersymmetric Higgs mixing 
mass is constrained by \9f^\ < vr/lO. 
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Figure Captions 



Fig. 1) Contours of constant ^^h!^ = 0.25, 0.5, and 1.0, as a function of mo and M, 
for Ao = 300GeV,^A = 0.8 tt,^^ = 0, and tan/3 = 2.1. The vertical lines 
represent contours of constant bino purity, p — 0.95 and p — 0.99. The 
shaded regions yield an LSP which is either a chargino or a sfermion. 

Fig. 2) Contours of constant Mmin, the minimum gaugino mass parameter needed 
to bring the neutron electric dipole moment below experimental bounds, for 
Aq = 300 GcV, mo = 100 CeV. The light central region corresponds to Mmin < 
200 CeV, and successive contours represent steps of 100 CeV. Note that M < 
400 CeV implies 9^ < n/15. 

Fig. 3) Same as Fig. 2 for the electron electric dipole moment. In this case M < 
400 CeV implies 9^ < tt/IO. 
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Figure 3 
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